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1. Introduction 

The study of quantum particles confined to manifolds of a mixed dimensionality has a 
long history starting from [TT] . Recently this problem attracted a new interest connected 
with possible influence of external fields and/or internal degrees of freedom. In this 
paper we continue this line of research and study the dynamics of a charged particle with 
spin i, having in mind an electron, which moves on a "hybrid surface" consisting of a 
halfline attached by its endpoint to a plane. In the plane the electron interacts with a 
constant magnetic field with orientation perpendicular to it, and with its own spin via 
a spin-orbit interaction — in the present work we suppose that the latter is of the form 
due to Rashba. 

The coupling at the contact point between the halfline and the plane can be chosen 
in different ways to make the resulting Hamiltonian self-adjoint. In this sense it may 
include a point interaction; such a motion in the plane alone was analyzed in [3J; the 
hybrid plane without the magnetic field was discussed |13j . 

The construction of the Hamiltonian follows the usual pattern using the theory of 
self-adjoint extensions. As a starting information we use free evolution on the discon- 
nected parts of the configuration space. On the one hand it is the halfline equipped 
with the Laplacian with Neumann boundary condition. This is, of course, simple. More 
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complicated is the other part describing a particle with spin moving in the plane under 
the effect of a constant magnetic field and Rashba interaction — here we can use the 
results of the papers [5] and [S] which we briefly recall in the next section. 

After this preliminary we construct in Section 3 Hamiltonians describing the coupled 
system; we write the appropriate generalized boundary conditions and derive the corre- 
sponding Green function. In Section 4 we study properties of these Hamiltonians for a 
spin-independent coupling. We analyze scattering of a particle travelling on the halninc 
finding the corresponding reflection amplitude. Furthermore, we find the spectrum of 
the Hamiltonian as well as the resonances the origin of which are the Landau levels in 
the decoupled plane. 



2. A preliminary: motion in the plane 

In this section we collect some results about the motion in the plane needed in the 
following; proofs and details can be found in [5]. We consider a charged two-dimensional 
particle 1 of spin \ under the influence of a uniform magnetic field of intensity B orthog- 
onal to the plane — for definiteness we suppose that B > 0. Let A be the corresponding 
magnetic vector potential, B = -g^- — ^^; we employ the symmetric gauge putting 
A(x,y) — (\By, ~\Bx). The Hilbert state space is 

^plane = L 2 (R 2 ,C 2 ) 

If we take into account the Rashba spin-orbit interaction, the motion of the particle in 
the plane is described by the following Hamiltonian 

Hr = H + j^Ur + ^ \i b Bg z (2.1) 

where 

1 , 

• Hq = II do with Hf := —ihdj — ~ Aj, j = x, y, where erg is the 2x2 identity 

matrix, 

• Ur — UxILy — (J y Ii x , with the standard notation for the Pauli matrices, a x = 

1 \ / -i A / 1 

1 )'°V = l i J'^= -1 



and hb = — is the Bohr magneton. Furthermore, m e and m* are the electron mass 
2m e c 

and its effective mass, respectively, g* is the effective g-factor, and ctR is the real-valued 
Rashba constant. 

The values of the units are not important in the following hence we shall use mostly 

dimensionless coordinates introducing the following notation: the coupling xr := — * R , 

nr 

lr To hedge us against a terminological objection we hasten to add that we have in mind a two- 
dimensional model of a real-world particle which proved to be extremely efficient, in particular, as a 
description of a rare electron gas confined to a very thin layer. 
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the magnetic flux quantum <J>n := , furthermore, b : = — B and a := — A or 

e $o $o 

a = (\by, — \bx). Setting now 7 := — ^ — we can rewrite the above Hamiltonian as 

H R = K 2 <t + 2h r U r + 7 &CT Z 

where K := ^ p — a and {/r = a x K y — a y K x . 

As usual the properties of such a Hamiltonian are encoded in its resolvent. The latter 
is known explicitly and setting 

PR ■= °> VR 



2xr 



^z + xl + Pl , C| (*>) := (VR ± + b - 0_ 



it is possible to write an explicit formula for the Green function components 

G£(x,x / ;*) Gi|(x,x';z) 



Gfl(x,x';z) = 
where x= (x,y), x' = (x',y'), and 



G£(x,x';z) Gf(x,x';z) 



G#(x,x';z) 



Pr- xr, 



-**(* v':-: 



G 2 fl 2 (x,x';^) 



(G (x,x';Q(b))-G (x,x';Ci(6) 
+ 2(G (x,x';C fi (fo))+G (x,x';Ci(fo))) 
'Go(x,x';Q(-6)) -G (x,x';CA(-&) 

+ i (G (x, x'; Ci(-6)) + G (x, x'; £ (-6))) , 
while the off-diagonal elements are 

= |6| ((x x>) - i(y y')) Q{ _ b)) 
-G (x,x';CA (-&))! + k(x,x';Cfl(-&)) - F (x,x'; 



and G|j 1 (x, x'; z) = G fl 2 (x', x; z). Here we have used the notation 



P o(x,x';z) * ( * - JW- 



4tt \2\b\ 



2 2\b\' 



xex P (f(x A xO-M (x _ x02 )^_^, 2; M (x _ xr) . 



|b| 



2|6| 



and 



xexpg(xAx')-^(x-x') 2 )^ 



2|6| ' ' 2 1 ; /' 



(2.2) 



(2.3) 
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where <J> is the confluent hypergeometric function and T is the Euler gamma function. 
As usual in these situations it is necessary to know also the renormalized Green function 
with the diagonal singularity removed given by 

G^(z)~ lim [G K (x,x';z)-S(x,x') 
where S(x, x'; z) := — log |x — x'|o- - If we put 
£/(,):= Jim (G (x,x';,)+i-log|x-x'|) =-i-(v-(i-^)-2V>(l)+logM), (2 . 4 ) 



then 

G5?"(z) 



G i? n i( z ) o 
G^ 2 (z) 

/ U(Q(b))-U(&)) 
2^ I U(( R (-b))-U(C+(-b)) 

W U(C R (b)) + U(C+(b)) 



M U(Q(-b)) + U((+(-b)) J > 



where 



_ J_ / 1 _ fe-x r \ / ^-(y/z + ^ + ^ + x r ) 2 \ (2 5) 



and 



R ^ 2 ^ ' ' Att \ 9 i /_ , Q2 i 772 \ ^ 



4n\2 2y/z + p 2 R + x? J \ 26 

i(l fa + *r \ / 26 + ^-( v /z + ^ + x2 + ^) 2 ^ 

4^2 fyz + ft + MtjVy ~ 26 

f 

4tt 



(2.6) 



The symbol ip(z) here means the digamma function, known to be a meromorphic function 
of z with no branch cut discontinuities and with simple poles at z = 0, — 1, —2, defined 
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by 

1 

%b(z) = -7 + (z - 1) ^ 7 T7 r . 

y ' v ' ^ (n + l)(z + n) 

n=0 y Jy ' 

The spectrum of the magnetic Rashba Hamiltonian has been derived in [5] . It consists 
of infinitely degenerate eigenvalues which are natural to call modified Landau levels, 

a p (H R ) = {e ± (n, s) : n e N, s = ±1} 

/ (2-7) 

e ± (n,s) = |6|(2n + l - s sign b) ± 2x R ^J (3 2 R + |6|(2n + l - s sign b) ; 

one can recover this result by inspecting the resolvent singularities in ( |2.3[ ). 

3. Motion in the hybrid plane 

After this preliminary we come to the problem described in the introduction and 
suppose that the configuration space consists of a plane described above to which a 
halfline lead is attached; without loss of generality we place the junction to the origin of 
coordinates in the plane. The construction is analogous to the non-magnetic case |13j 
the difference being in the plane Hamiltonian component, nevertheless, we describe it in 
sufficient detail in order to make the present paper self-contained. 



3.1. The halfline-plane coupling 

The lead Hilbert space is %i ca d = L 2 (R + , C 2 ), and the whole state space of the system 
is the consequently the orthogonal sum % := "Hicad ® 'Hpianc- I n other words, the wave 
functions are of the form '5 = {V'lcad, "0pianc} T with each of the components being a 2 x 1 
column. The construction starts from the decoupled operator H° :— i?i ca( j ® Hr where 
the first component is the Laplacian on the halfline iJioadV'ioad = — ^[cad with Neumann 
boundary condition at the endpoint 1 , while Hr is the magnetic Rashba Hamiltonian 
discussed in the previous section. Following the method of [TT] we restrict the operator 
H° to functions which vanish in the vicinity of the junction obtaining thus a symmetric 
operator of deficiency indices (4,4). In the second step we construct its self-adjoint 
extensions which are regarded as admissible Hamiltonians. 

The extensions can be characterized in various ways. There is a general scheme com- 
ing from von Neumann theory, however, it is by far more practical to employ boundary 
conditions when dealing with a problem of the present type. The boundary values on 
the halfline are the simply columns Viead(0+) and ip[ eSu< i(0+) ■ On the other hand, in 
the plane we have to use generalized ones. The functions in the domain of the restric- 
tion have a logarithmic singularity at the origin and the generalized boundary values 
Lj (V'piane), j=0,I, appear as coefficients in the corresponding expansion, 

V'planc(x) = -— L {lppl ane ) ln|x| +£l(V>pl aI1 e) +o(|x|). (3.1) 



1 There is no need to consider the effect of the magnetic field on the halfline even if the lead is not 
perpendicular to the plane because one can always remove it by a simple gauge transformation. 



[Author and title] 



6 



where 

io(V'plano) = lim - 

|x|->0 ln|xj 

il(^plane) = lim [V>plane(|x|) - ln(|x|)L (^plane(|x|))] 
|x|->0 

Using these boundary values we can write the sought boundary conditions as 

Cad(0+) = ^cad(0+) + C*L O (^planc), ^ ^ 

il(V'planc) = Clp lcad (0+) + DL (t/j plaDC ) , 

where A, D, C are 2x2 matrices, the first two of them Hermitian, so the matrix A := 
{ J q C £ ) ) depends of sixteen real parameters as the deficiency indices suggest. One can 



check easily that the corresponding boundary form vanishes under the condition (3.2), 
which means that each fixed A gives rise to a self-adjoint extension of the restricted 
operator. 

It is worth noting that the above boundary conditions are generic but do not cover all 
the extensions leaving out cases when the matrix A is singular; this flaw can be mended 



in the standard way |17|,f3] if one replaces (3.2 ) by the symmetrized form of the relation, 



</w(o+)\ , ^V4ad(o- 



V-kuCVWue)/ + V-MVVanc)/ 

where A,B are matrices such that (-4|£>) has rank four and AB* is Hermitean. We will 



restrict ourselves, however, to the generic case B = —I expressed by (3.2 1 in the following; 
the same is true for the alternative form of the b.c. mentioned below. 

The way in which the parameter matrix is chosen depends on physical properties of 
the coupling between the lead and the plane. In particular, diagonal A, D, C correspond 
to the situation when the junction does not couple the spin states, and moreover, scalar 
matrices describe a spin-independent coupling. It is obvious that the lead and the plane 
are decoupled if A is block-diagonal, i.e. C = 0. A naive interpretation of the conditions 



(3.2 1 is that C is responsible for the coupling while A and D are point perturbations at 



the two components of the configuration space, respectively. 

Remark 1. An attentive reader may wonder about the relation between the boundary 
conditions and the geometry of the problem, that is, the angle between the halflinc 
and the plane. Likewise, one may ask whether the coupling could be influenced by the 
magnetic field. The above analysis gives no answer to these questions; it only guarantees 
self-adjointness of the Hamiltonian, or in other words, conservation of the probability 
current through the junction. A natural approach would be to consider a "fat" hybrid 
plane in analogy with the analogous problem for networks - cf. |10j and references therein 
- and to analyze the limiting behavior as its thickness tends to zero. This question is 
open and by far not easy. 
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3.2. The Green function 

As usual properties of an operator are encoded in its resolvent, hence our next task is 
to find the latter for the above constructed self-adjoint extensions. A suitable tool to do 
that is Krein's formula [21 U\ which allows us to find the sought resolvent starting from 
Green's function of the decoupled system which has a block-diagonal form, 

[x,x ,x,x ,z) — y Q2 G R (x,x';z) ) ' {SA) 



where O2 is the 2x2 null matrix, G,r(x, x';z) is given by (2.3 1 and 



G\ ca d(x, x';z) = —j= cos y /zx < e l ^' x> a 
V z 

with the conventional notation, x < := min{a;,x'}, x> := max{a;, x'}, since we have 
assumed Neumann boundary condition at the halfline endpoint. We introduce Krein's 
function Q(z) , which is an analytic 4 x 4-matrix valued function of the spectral parameter 
z, as the diagonal values of the kernel, with the above described renormalization in the 
planar part, specifically 

9(*):=( J^V (3-5) 



2 G%«(z) 

The full Green function is obtained by a finite-rank perturbation of the free one. It is 



convenient to rewrite the conditions (3.2 1 using a modified basis in the boundary value 



space: instead of the vectors employed above we take 

f ^ s= (-^(o+A f ^ := ^w(o+) 



an 

One can check easily that they satisfy ATitp + BT21P — with B — —I and 

1 - ( -a- i c* \ 

A.-y _ CA -i D - C A- l C* ) ■ 

It is clear that A which equal up to the sign to AB* is Hermitean. The reason for the 
modification is that with the last boundary conditions our comparison operator H° is 
characterized by Tiip — 0, i.e. A = I, B° = 0. This allows us to use the result of [3] 
directly (in our case it is nothing else than the usual Krein's formula) to infer that the 
resolvent kernel of is given by 

Ga(x,x';x,x';z) = Gr{x t x'; x, x'; z) (3.7) 
-G°(x, 0; x, 0; z) [Q(z) - A] _1 G°(0, x';0, x'; z) ; 

the second term which is a rank sixteen operator represents the resolvent difference 
between the free and full Hamiltonian. It is important to note that even in the situation 
when the coupling is spin-independent, A = (^) ®ffo and similarly for A, the Green 
function does not decompose; the reason is that the spin states are still coupled by the 
spin-orbit interaction in the plane. 
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4. Properties of 



In this section we will investigate the coupling between the halfline and the plane 
with a particular choice of the boundary conditions containing a parameter which allows 
us to control the coupling strength. Specifically, we employ the following choice of the 
matrices, 



A e = 



a e 
e d 



' cr ■ 



e^O and B £ = - 



1 
1 



(4.1) 



so the two parts of the configuration manifold are coupled and there is no coupling 
between spin degrees of freedom at the contact point. In the plane, of course, we have 
the spin orbit interaction. In the decoupled case, e — 0, the constants a, d can be regarded 
as point interaction strengths at the halfline endpoint and in the plane, respectively, to 
which a physical interpretation can be given in analogy with |12j. The off-diagonal terms 
in the matrix A e are responsible for the coupling between the halfline and the plane. 

Having chosen the coupling our next task is to specify the quantities appearing in the 
general formula (3.7); we get 



(Q(z)-AT 1 



De{z) 



( 



r§,i(*) 
V o 



ri, 2 (z) 
o 

ri, 2 (z) 



rf, a (*) 



r§, a (*) 







ri, 4 (z) 
o 

Oi(*0 J 



with 



D e (z 

rf,i(* 
rf >3 (* 
r| >a (* 



(Q e "W- 



= n 2 (z) = £ 



\7z 



i 

Tz 



V~z 

G^; 2 (z)-d 

G r &(z)-d] -& 



G^{z)-d 



— e 



Gra(z) - d 



— e 



V~z 



V~z 



G x i%(z)-d 
G%%(z)-d 
G^(z)~d 



— e 



— e 



4.1. Scattering 

The analysis of transport on a hybrid surface of mixed dimensionality follows the 
scheme described in numerous papers — see, e.g., [TT], [13], [T5], [5], [5], and references 
therein. 
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If we have a scattering system in which the difference between the two dynamics 
is expressible in terms of self-adjoint extensions of a symmetric operator with finite 
deficiency indices it is possible to write the scattering matrix directly in terms of the 
Krein operator- valued function Q(z) and the boundary conditions - cf.pQ or [3]. Using 
this approach — we refer specifically to formula (1.5) from [3] — we arrive at 



/ (-7?- a )( G fi°i( z )- d > £2 
(^- 5 )( G « on i( z )- d > £2 











-a)(G r H ep 2 (^)-rf)- 



V 



— r: 

-a) (G™° (z)-d)-e* 

1 

1 / 



(4.2) 



Putting z = k 2 we can in analogy with [T5] describe the reflection amplitude of a particle 
travelling along the halfline, with the "upper" spin component only, as follows 



K(k) 



5) (GfiCfc 2 ) 



d )-e 2 



a) [G^(k 2 



(4.3) 



which looks like the result contained in 13J. As briefly mentioned in the said paper, 
however, the magnetic case differs substantially from the non-magnetic one. The point is 
that the Green functions G^"(fc 2 ) with j = 1,2, are real-valued, and as a consequence, 
the scattering on the halfline is unitary, |7?.(fc)| 2 = 1, in the magnetic case. From the 
spectral point of view one can expect formation of resonances due to the perturbation 
of the discrete spectrum of the spin-orbit Hamiltonian in the plane embedded in the 
continuous spectrum of the free Hamiltonian on the half-line as it is common in similar 
cases — see again [14], [15] and references therein. 

4.2. Spectral properties of iJ- 4 

The literature devoted to the effect of weak or local perturbations on the Landau 
Hamiltonian is rich. In particular, the recent paper 48] investigated a weak perturbation 
of the Landau Hamiltonian by a fast decaying or even compactly supported electric 
or/and magnetic field; the effect on the Landau levels was gcncrically a splitting. Of 
course, eigenvalues which split off a Landau level appear in many situations — recall, 
e.g., the classical analysis [16] of a point interaction in a homogeneous magnetic field, 
possibly in presence of a potential. However, the result of [18] was more surprising: 
the authors discovered that a 2D axially symmetric short-range potential gives rise to 
an infinite number of the negative-energy levels if one takes into account the spin-orbit 
interaction. 

Spectral properties of require an accurate analysis, in particular, because the 
decoupled system has numerous embedded eigenvalues, and while Weyl's theorem guar- 
antees stability of the essential spectrum, its character might change by the perturbation. 
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For the sake of definitcness we suppose the magnetic field intensity is positive, b > 0, in 
the chosen coordinate frame. 



Theorem. Assume that a > 0; then the spectrum of Hj^ looks as follows: 

(a) the point spectrum o~ p (H^c) = o~ p (Hr) U E 7 where £ is a finite set of negative 
eigenvalues, 

(b) the continuous spectrum a c (H^) = [0, oo), 

(c) there are infinitely many resonances with the real parts in the gaps between the 
eigenvalues of Hr and negative imaginary parts. If the coupling is weak, their dis- 
tance from the embedded eigenvalues corresponding e — is 0(e 2 ); the perturbation 
expansion is given by the relations (4- 14) below. 

Proof. The Weyl theorem guarantees the preservation of the essential spectrum, 
o~ ess (HA e ) = a ess (Ho). Since both the operators have a common symmetric restriction 
with deficiency indices (4,4), it follows from general principles [19j Sec. 8.3] that the 
negative spectrum of A e consists of at most sixteen eigenvalues, multiplicity taken into 
account. To learn more about the the point spectrum, o~ p (H^), it is necessary to check 
directly the singularities of the resolvent (3.7 1. We will analyze explicitly one component, 
the other ones can be treated similarly. We rewrite the (3, 3) component of the matrix 
(E7| as 



(G A (x,x' ';x,x' '; z))( 3 ,3) 



Grt( x i x 'i z ) 



G 



£(x,0;*)(^-o) G%(0,x?;z) 





Cfl( X 7 X 'i Z ) 


(*-*) 


G*%%z)-i 


-e 2 



where G}j(x, x'; z) equals 

1 Pr - 

2 2rj R 



G (x,x';Cfl(&)) 



1 _ Pr - xr 

2 2 m 



G (x,x';Ci(&)) 



(4.4) 



Expanding the digamma function in Go (x, x'; z) we get 



G (x,x';z) = ^e(x,x') 



r=0 



2b) r 



1 2 2|6| ' ' 2 K ' 



l ( ,,(M (x _ x > 



-4tt U(z - 2|6|r) - log ^ + 2*(1) - 2*(1 + r) 
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where 



and 



9(x,x') :=cxp(|(xAx')-^(x-x 



<I>(a, c, x) 



(a) r 



oo 

S ( c ) rH 

r— 



x r with (z) r 



T(z + r) 
T(z) 



It is easy to check that the first term of the last expression at the rhs of (4.4 1 is still 
singular for Cb(&) = n an d these singularities are the same as those of the decoupled 



Hamiltonian. In a similar way we can treat the other components of the resolvent (3.7| 



concluding that the Rashba energy levels remain to be embedded eigenvalues of infinite 
multiplicity in the coupled case. Using this fact and the above mentioned preservation 
of the essential spectrum it is possible to determine the continuous part of the spectrum, 



The resolvent (3.7) can have naturally other singularities coming from zeroes of D e (z). 
In the decoupled case, e = 0, the point interaction on the halfline has one bound state 
with the negative energy Eq = —-gi, since we have assumed a > 0. Let us analyze the 
effect of the coupling, e^O, observing the component of (3.7) with j = 1, 2, 



{Ga e (x, x ; x, x'; 



cos V zx<- e 



-iWzx ~—l\/zx 



The only singularity in the vicinity of Eq is obtained from the equation 



J = l,2 



The Rashba eigenvalues of the unperturbed Hamiltonian, corresponding to s = 0, d = 0, 
are all real, positive and infinitely degenerate, and since the perturbation is of finite rank 
in the resolvent sense, continuity of the singularities w.r.t. e is guaranteed. 

For a start let us discuss qualitatively the solutions of the equation D e (z) — around 
Rashba eigenvalues. For the decoupled case, e — 0, and d ^ we a point interaction 
in the plane with spin-orbit interaction and the effect of this perturbation is analogous 
to that of pR>| . Suppose that d is sufficiently large to keep all the eigenvalues real and 
positive — we recall that the lowest one goes to — oo as d — > — oo — which means in 
the language of two-dimensional point interactions that the perturbation is sufficiently 
weak. Switching then the coupling on, e =^ 0, the condition D £ (z) — with the solution 
z = X G R+ can be then rewritten as 



Vx 



G33(A)-d 



GS3(A)-d 



1.2 



(4.5) 
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and we notice that this equation cannot be solved by any real positive A because on the 
left-hand side we have always a real positive number while the right-hand side is always 
purely imaginary. This is due to the fact that the function G^"(A) takes real values for 
every A > and j = 1,2. If there is a solution to this equation it must be therefore 
outside the real axis. 

On the other hand, for a sufficiently strong point interaction in the couple case, e = 0, 
in other words for d sufficiently large negative there may be negative eigenvalues. The 
equation for singularities in this case can be written as 

G%(z)=d (4.6) 

From the properties of the digamma function we see that the function G^"(z) is 
monotonously decreasing around the origin, it is positive in a neigborhood when it ap- 
proaches the first Rashba level and has a zero on the negative halfline. Consequently, 
for d sufficiently large negative there is a negative solution. Switching now again the 



coupling in, e ^ 0, an taking z — —A with A > we see from (4.5) that such a solution 
is still a negative eigenvalue, at least for |e| small enough. 

Let us now look more closely how the singularities of the resolvent behave in the weak 
coupling regime. We are going to find the corresponding series expansion in e using a 
recursive procedure. To begin with we rewrite the equation D £ {z) = as 

1 

where we have introduced 



/(*) = 7^73 (4-7) 



f( z ). = 1 = {ad - e 2 )ay/z + d ie 2 ^fz 

d- {ad-e 2 ) 2 z + d 2 {M-e 2 ) 2 z + d 2 

To simplify the recursive procedure we adopt the following notation, 



Ef n = 2nb - 2k t a//3| + 2n6 
E$ n = 2nb + 2>t r ^l3 2 R + 2nb 



E° n = 2{n + 1)6 - 2x r ^/p 2 R + 2{n+ 1)6 
E% n = 2{n + 1)6 + 2x rV /(3 2 R + 2{n+l)b 



(4.9) 



and we will look for the fixed points of 

^^(BMfzW-^)^/^- 1 )), 1 = 1,2, 3, 4, j = 1,2, (4.10) 



G^-(zW) 

where 



P(^n)=«T T^TTld . 1 = 1,2. (4-11) 



G^{z) 



=E» 
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PiKn) = 



dz 



I = 3,4. 



p(4») = ^ 3 Vi) = - 

p(E? n ) = = - 



8tt 



b p 2 R + y/2nb + [3 
8tt ^2nb + 



(4.12) 



(4.13) 



In this way we can determine the position of the singularities in the leading order in e. 
Starting with z^ = Ef fl + p ^ ( , > j and / = 1,3 we have two negative real solutions if 

d^mm{C 1 ,C 2 } with d = (E^piE^))- 1 and C 2 = (^iPCEg,!)) -1 . 

On the other hand, for d sufficiently large the starting point of the recursive procedure 
will be a positive eigenvalue z^ ' = E® n + p ^ ^ i and the fixed points of the recursive 



procedure, from (4.10) and (4.13), will have a positive real part and a negative imaginary 



part; in leading order in e the real and imaginary parts of the resonances are 

1 



TZ(E t 



■res \ 
Ln ) 




(4.14) 



p(EZ n )d*(l + a*E° ) 



From the first of (4.14), putting e = 0, it is obtained the displacement of the Rashba 



energy levels due to the point interaction in the plane. In the coupling case, for e 7^ 
the energy eigenvalue will migrate in a resonance with negative imaginary part. 
The convergence of the sequence in a ball of radius e 2 centered at z^ = E? „ +- 

can be proved using the following estimates, 



,(1) 



,(°)ll - 



p(E? n : 



ll/(- (0) )ll < 



|| z (fc+i)_ z (fe)|| ^c 2 s 2 \\z^-z^\\, 



where the constants Cj are independent of s. 



□ 



Remark 2. The assumption a > was used to ensure that there are proper eigenvalues 
in the decoupled case on the halfiine. In a similar way one can treat the case a < 
where for e = there is a singularity for the decoupled resolvent but it is on the second 
Riemann sheet, i.e. an antibound state. One can check that for e ^ this point will 
remain on the real negative axis on the second Riemann sheet. 



In order to illustrate the above conclusions, one can analyze in an example the effect 
of the coupling on the positive eigenvalues coming from the Rashba energy levels in the 
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plane. For instance, let us consider the following eigenvalue in the decoupled case: 

E° ia = E 3A = 26 - 2x RV //3| + 2& (4.15) 

We fix the value of d = 1 and study numerically the effect of the coupling e 7^ 0. In the 
following pictures the real and the imaginary parts of the resonances arising from the 
same Rashba eigenvalue 2 = E® 1 are shown. 




Figure 1: In blue the real part of the resonance is plotted as a function of e, coming from 
El 2 (color online). In red (the lower curve) the real part of the resonance coming from 
Eg J with d = 1, 6=1, 5 = and x/6 = 0.1 is shown. 
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Figure 2: The imaginary part of the same resonances as above. 

As it often happens in such resonance problems — see, e.g., [T3] — for large values 
of the coupling constant e the resonance can approach an eigenvalue again. In our case 
it is clear that |3.2| reduces in the limit e — > oo to Dirichlct boundary conditions. 
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